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ABSTRACT 

Constraints on the clustered mass density fi m of the universe derived from the 
observed population mean intracluster gas fraction (//cm) of X-ray clusters may be 
biased by reliance on a single-phase assumption for the thermodynamic structure of the 
intracluster medium (ICM). We propose a descriptive model for multiphase structure 
in which a spherically symmetric ICM contains isobaric density perturbations with 
a radially dependent variance <J 2 (r) = 0^(1 + r 2 /r 2 )~ e . The model extends the work 
of Gunn & Thomas (1996) which assumed radially independent density fluctuations 
thoughout the ICM. 

Fixing the X-ray emission profile and emission weighted temperature, we ex- 
plore two independently observable signatures of the model in the {er c , e} space. For 
bremsstrahlung dominated emission, the central Sunyaev-Zeldovich (SZ) decrement in 
the multiphase case is increased over the single-phase case and multiphase X-ray spec- 
tra in the range 0.1 — 20 keV are flatter in the continuum and exhibit stronger low 
energy emission lines than their single-phase counterpart. We quantify these effects for 
a fiducial 10 8 K cluster and demonstrate how the combination of SZ and X-ray spec- 
troscopy can be used to identify a preferred location {er c , e} in the model plane. From 
these parameters, the correct value of (//cm) in the multiphase model results, allowing 
an unbiased estimate of fi m to be recovered. 

The consistency of recent determinations of the Hubble constant from SZ and X-ray 
observations with values determined by other methods suggests that biases in ICM gas 
fractions are small, < 20%. 
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1 INTRODUCTION 

The existence of non-baryonic or "dark" matter on very 
large scales in the universe is inferred from a number of 
observations, including X-ray and gravitational lensing ob- 
servations of galaxy clusters. Observations suggest that the 
baryonic component of clusters predominantly consists of 
hot, diffuse, intergalactic medium (ICM) which emits X-rays 
by scattering of electrons in the Coulomb fields of electrons 
and ions, i.e., thermal bremsstrahlung. The X-ray observa- 
tions determine the ICM mass content in a model depen- 
dent fashion. Recent analysis of the flux limited Edge sam- 
ple employs the standard, isothermal /3-model and finds a 



mean ICM mass fraction {fichi) = 0.212 ± 0.006 (Mohr, 
Mathiesen & Evrard 1998; see also White & Fabian 1996; 
David, Forman & Jones 1996) within the virial regions of 
27 nearby clusters with X-ray temperatures above 5 keV. 
This value is several times larger than that expected in an 
Einstein-deSitter universe with the observed light element 
abundances (White et al. 1993). 

One way to reconcile the cluster observations with a 
universe having critical mass density Q m = 1 is to suspect 
that the standard model treatment of the ICM posseses sub- 
stantial systematic errors. Gunn & Thomas (1996, hereafter 
GT96), motivated by models of cooling flows (Nulsen 1986; 
Thomas 1988) propose that a multi-phase ICM structure ex- 
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ists throughout the cluster atmosphere. A given macroscopic 
volume element contains gas at a range of densities and tem- 
peratures which are assumed to be in pressure equilibrium. 
Fixing the gas mass within this volume, the emission mea- 
sure of a multiphase gas will increase as the clumping factor 
C = (p 2X ]/(p^ . But since we observe luminosity, not gas 
mass, the implication is that clumped gas requires less total 
mass M gas oc 1/yC in a given volume to produce a fixed 
X-ray emissivity. 

The standard analysis of the cluster plasma assumes 
that it exists in a single thermodynamic phase at any loca- 
tion within the cluster. In most cases an isothermal, "beta" 
model (Cavaliere & Fusco-Femiano 1976) is used to describe 
the cluster plasma electron density for a spherically sym- 
metric atmosphere. Under these assumptions, the observed 
azimuthal X-ray surface brightness profile determines the 
volume emissivity at radius r from the cluster center 



^r) = p 2 (r)A x (T x )=^[l+ r - 



-3/3+4 



(1) 



Here £o is the central value of the X-ray emissivity, r c is 
the core radius of the X-ray emission, Ax(Tx) is the (suit- 
ably normalized) plasma emission function at a tempera- 
ture Tx over a prescribed X-ray bandwidth. The temper- 
ature Tx is determined from X-ray spectral measurements 
by, for example, fitting the observed spectrum to a thermal 
bremsstrahlung model. With observations and plasma emis- 
sion model in hand, one then constructs the gas mass den- 
sity p(r) = (£o/Ax(T x )) 1/2 (l + r 2 /rl)- 30/2 and integrates 
outward from the origin to define enclosed gas mass. 

The total (baryonic plus non-baryonic) mass within a 
radius r is inferred from assuming that the plasma is in hy- 
drostatic equilibrium, supported against gravity entirely by 
thermal pressure. The fluid equation of hydrostatic equilib- 
rium then sets the total, gravitating mass 

M "W = -GpW> (2) 
which for the fiducial, single-phase, isothermal /3-model 
cluster gives 



M toM (r) 



3/3 ksTx 



Vr 



G pm p l + r 2 /r| 



(3) 



In this paper, we extend a multiphase ICM model first 
proposed by Gunn & Thomas (1996, herafter GT96) to in- 
corporate radial variability in the multiphase structure. Ra- 
dial variability is a natural expectation. Since both cooling 
timescales (oc p _1 if nearly isothermal) and local gravita- 
tional timescales (oc p -1 ^ 2 ) increase outward from the cluster 
core, the timescale for development of multiphase structure 
should also be larger at the virial surface than in the core of 
a cluster. 

We introduce the theoretical model in §2 below. In §3, 
we examine the effects of a multiphase structure on the 
mean intracluster gas fraction (fioM) inferred from X-ray 
observations and consider observable implications for the 
Sunyaev-Zeldovich (SZ) effect and X-ray spectroscopy of 
the ICM. For the latter, we examine two specific signatures 
— the excess (relative to the single-phase case) in central SZ 



decrement and an X-ray spectral hardness ratio — for the 
case of a "Coma-like" cluster. We show how the pairing of 
X-ray spectroscopy and SZ image can be used to estimate 
the magnitude of systematic error introduced into estimates 
of ICM gas fraction by assuming the standard /3-model . 



2 THEORY 

GT96 argue that if a spectrum of plasma density fluctua- 
tions were generated in a cluster at a substantial fraction of 
a Hubble time in the past, then its densest phases would cool 
and be removed from the plasma. Following Nulsen (1986), 
they argued that this would produce a power-law spectrum 
of fluctuations f(p) oc p -7 at the present, formed from a nar- 
row range of phases that were initially tuned to have cooling 
times comparable to a Hubble time. However, this argument 
ignores the stochastic nature of gravitational clustering in 
hierarchical models of structure formation. In such models, 
clusters grow largely by mergers of proto-cluster candidates 
embedded within the large-scale filamentary network. It is 
suspected that strong mergers may, through plasma tur- 
bulence, effectively "reinitialize" density fluctuations in the 
ICM. Since the time since the last major merger is a random 
variable in a coeval population, then a volume limited sam- 
ple will contain clusters whose multi-phase structures are at 
different stages of development. This idea is consistent with 
observed properties of the local X-ray cluster population, 
in which a range of central cooling flow behavior is present 
(Fabian 1994). 

Because of this and other uncertainties in the dynamical 
development of multiphase structure, we postulate a log- 
normal form for the multiphase density perturbations. We 
do not attempt a formal justification for this choice; it is 
motivated largely by a condition of "reasonableness" and 
the fact that it simplifies calculations below. The formalism 
requires only low order moments of the distribution, so the 
model can be recalculated for arbirtary f(p). 

We postulate the existence of plasma density fluctua- 
tions in a spherically symmetric cluster atmosphere which: 
(i) are isobaric at a given radius, (ii) produce a volume emis- 
sion profile consistent with equation ([[]) and (iii) exhibit 
an isothermal emission weighted temperature with radius. 
The first item is based on a hydrostatic assumption and the 
remainder impose observed constraints on the X-ray im- 
age and emission weighted temperature profile. Although 
isothermality extending to r2oo — the radius within which 
the mean total mass density is 200 times the critical density 
— may not be supported by observations (Markevitch et al. 
1998) or simulations (Frenk et al. 1998), temperature drops 
of only 10 — 20% are allowed within r2oo/3 (Irwin, Breg- 
man & Evrard 1998). Since the observables we stress in the 
analysis are core dominated, our results are not particularly 
sensistive to departures from isothermality which may exist 
near r2oo- 



Multiphase Intracluster Medium 



2.1 The multiphase distribution function 

We assume a log-normal form for the cluster plasma density 
phase distribution f(p) dp, the fraction of a volume element 
at a radius r that contains plasma of density of between p 
and p + dp, 



f(p)dp 



1 



■ exp 



In 2 [p/ P0 (r)] \ dp 
2a 2 (r) J p ' 



(4) 



The quantity po(r) is a reference density and a 2 (r) is the 
variance of the distribution. Since the core radius presents 
a characteristic scale in the X-ray image, we take a form 



(5) 



for the variance, with r c the core radius of the beta-model 
density profile described earlier, and a c and e are free pa- 
rameters which set the magnitude and radial dependence of 
the multiphase structure. We consider such a parameteri- 
zation in order to couple the magnitude of density fluctua- 
tions to the likelihood that the local conditions have allowed 
cooling to amplify them. A simple parameterization is one 
in which the variance scales with the inverse of the local 
cooling time cr 2 (r) ocr^^r). An isothermal atmosphere (for 
which T coo i(r) oc p -1 ^)) will have cr 2 (r) oc p(r), implying 
e = 3/3/2 = 1 for the characteristic (3 = 2/3 value seen in 
X-ray images. We consider values in the range e £ — 1. 
The limit e — > oo represents a multiphase structure exist- 
ing purely within the cluster core. In the limit a c — > 0, we 
recover a single-phase plasma for any value of e, while the 
limit e — > yields the multiphase model results (no position 
dependence) of GT96. 

The definition of po(r) is now absorbed into the speci- 
fication of the mean density at radius r 



(6) 



(p( r )) = I Pf(p)dp = p exp[ -a (r) 



where (} represents an ensemble average of volume elements 
on a spherical shell of radius r. A useful equation is a gen- 
eralization of equation to higher moments, namely 



(p") = / p q f{p)dp=(pY exp 



<?(<? 



1) 2, v 
cr (r) 



(7) 



2.2 A multiphase "isothermal (3— model " cluster 

We now impose some observational constraints on the 
model. Assuming a power-law emissivity function 



Ax(T)ocT Q , 



(8) 



the requirement that the emission weighted temperature 
profile be isothermal at temperature Tx implies that the 
condition 

(p 2 T 1+a ) 



Tx 



(9) 



<P 2 T Q > 
holds at all cluster radii. 

Under an ideal gas assumption P = (p/ pm p )kBT, with 
m v theproton mass and p the mean molecular weight, equa- 
tions (g) and (^) can be used to define the local gas pressure 
in the multiphase medium 



P(r) = *^(p(r)} exp[(l - a)a 2 (r)]. 



(10) 



pm p 

We now equate the known emission profile of the cluster, 
equation (|lj) , to the ensemble-averaged value of the emissiv- 
ity 

„2 



Aq / pUhp \ a . 



(11) 



The combination of equations (^|) and ( pd| ) is the canoni- 
cal "isothermal /3-model " assumption. From an observer's 
perspective, a multiphase cluster in these two measures is 
indistinguishable from the single-phase case. Equation ( |ll| ) 
can be rearranged to give 



(P(r)> 



C ™2 \ 1/2 

t,om p 
Ax(T x 

(a -!)(<* + 2) 2 



exp 



a 2 (r)). 



(12) 



Although this now defines the characteristic density po used 
in equation (^), it is better to identify the limit <r 2 (r) — > 
as the single-phase density. Following GT96, we introduce a 
multiphase "correction factor" for the gas mass C p (r) which 
relates the mean gas density in the multiphase case p m (r) 
to its single-phase value p s (r) 



p m (r) = (p(r)} ee Cp(r)p s (r). 
Equation ( Jl~2| ) then implies 

C p (r) = exp(^ ^ '-a (r)j. 

A similar exercise for the gas pressure 
P(r)ECp(r)(M£) Ps(r) 

yields 

C P (r) = exp( a (r)j. 



(13) 



(14) 



(15) 



(16) 



Note that for values of the X-ray emission exponent a < 1, 
the multiphase gas mass is lower than that of the single- 
phase model while the multiphase pressure is greater than 
the single-phase pressure. This arises because the high- 
density phases are more efficient in producing a given X-ray 
power (provided the emission is only a weak function of tem- 
perature). Since the emission weighted Tx reflects the tem- 
perature in higher than average density regions, the pressure 
at all radii is increased over the single-phase case. 

The cluster gas mass for the multiphase model within 
a radius r is given by 



M gas , m (r) =4vr / C p (r') p s (r') r U dr', 



(17) 



so that the enclosed gas mass for the multiphase model dif- 
fers from the single-phase case by the factor 

-^gas,m (r) 



C gaB (r) 



M gas , s (r) 

r C p {r') Pa {r')r' 2 dr' 
Jo Ps(r')r> 2 dr> 



(18) 
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Figure 1. The gas mass multiphase correction factor C ga s(0> 
equation , is plotted as a function of the dimensionless radius 
r / r 200 assuming parameters r c /r2oo = 0.1, a = 1/2 and j3 = 
2/3. The solid and dotted lines represent £ = 5 and e = 1.0, 
respectively, and values of cr c are as indicated. 



The total mass of the cluster within a radius r is deter- 
mined from the hydrostatic equilibrium, equation (^), when 
comibined with equations (|l2|)-(|l6|) give 



M to t,m(r) 



C ^ ( G P (r)P' s {r)+Ps{r)C' P {r) 
G 



Cp(t) 
C P (r) 



C P (r)p s (r) 
IcbTx 



G firrip 



C' P {r) 



C P (r) 



(19) 



The total cluster mass for the multiphase model differs from 
that of the single-phase model by the factor 



Ctot(r) 



+ 



k B T x 



Cp(r) 

Cp(r) ' GA'/ to t, s (r') pm p 



C' P (r 



C P (r) 



(20) 



For bremsstrahlung emission (a ~ 0.5), the gas mass is 
decreased and the total mass increased in the multiphase 
case, implying the enclosed gas fraction at radius r is lower 
than that for a single-phase medium by the factor Cb(r) = 

CgaeW/Ctotfr). 

Figure || plots the correction factor for the cluster gas 
mass, C gas , for a few choices of the controlling parameters 
(t c and e. For purposes of illustration, we take structural pa- 
rameters representative of rich clusters, namely r c /r2oo =0.1 
and /3 = 2/3 (Neumann & Arnaud 1999), and assume pure 
bremsstrahlung emission, a = 1/2. The effect of the radial 
falloff of the multiphase structure on gas mass estimates is 
substantial. Density variations with large central rms per- 
turbations a c ~ 2.0 produce substantial (factor ~3) relative 
correction to the gas mass in the cluster core, but the effect 
on the total virial gas mass (mass within r2oo ) is reduced 
to 25% if e = i and only 10% if e = 1. Degeneracies exist 
in the virial gas correction factor; a relatively weak mul- 
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Figure 2. The total mass multiphase correction factor Ctot(r), 
equation ( pp| ) is shown as a function of radius, using the same 



parameters and format as in Figure |l]. 



tiphase plasma distributed throughout the cluster can pro- 
duce an effect that is similar to a plasma with strong density 
variations concentrated toward the center of the cluster (cf. 
{<r c ,e} combinations of {1, t?} and {2, 1}). 

The correction factor for the total cluster mass, Ctot, 
for the same multiphase parameters is shown in Figure ^| 
By steepening pressure gradients, the multiphase effects in- 
crease the total cluster mass derived from equation (|l9[). 
Once again, weaker multiphase effects distributed through- 
out the cluster can yield a total mass within r2oo that is 
similar to a cluster plasma with stronger density variations 
concentrated in the cluster center. However, such concen- 
trated multiphase effects will produce a total mass profile 
that is steeper (c/.(<r c ,e) of (2, |) vs. (2,1)). Observations 
of strong gravitational lensing could be used to break this 
parameter degeneracy, to the extent that the hydrostatic 
assumption is valid in the cluster core. 



3 CONSEQUENCES 

We now turn to the issue of the effect of multiphase struc- 
ture on inferred ICM gas fractions and cluster observables. 
For the latter, we consider the effects of multiphase plasma 
on a cluster's X-ray spectrum and the Sunvaev-Zel'dovich 
microwave decrement through a line-of-sight taken through 
the center of the cluster. 

All of the results we discuss assume a standard structure 
model with core radius for the broadband X-ray emissivity 
(equation ([!])) of r c = 0.1r2oo, exponent f3 — 2/3, and, for 
creation of X-ray spectra, an emission-weighted X-ray tem- 
perature Tx = 10 s K. Unless otherwise stated, we employ 
a value of a = 0.36 for the exponent of the plasma emission 
function, derived from a Raymond-Smith code as described 
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Table 1. Model Definitions 



Label 


Cc 


e 


SP 







MP-A 


0.5 





MP-B 


2.0 


0.8 


MP-C 


1.0 





MP-D 


2.0 


0.4 


MP-E 


2.0 






below. Since we ignore galaxies in our modeling, the ICM 
gas fraction is synonymous with the cluster baryon frac- 
tion. We use the terms interchageably below, but it must 
be remembered that the stellar content of cluster galaxies 
and intracluster light presents an absolute lower limit to the 
baryon content of clusters. 

3.1 Baryon fraction bias 

The effects of increased total mass and decreased gas mass 
shown in Figures |l| and ^| combine multiplicatively to re- 
duce the cluster baryon fraction. The magnitude of the ef- 
fect within the virial radius is shown in Figure ^, where we 
show contours of Cb(r2oo), the baryon reduction factor, in 
the {er c , e} plane. 

The baryon reduction effect peaks at high a c and low 
e. At e = 0, the uniform, mulitphase structure of GT96 is 
recovered, with magnitude 



91 

C P 



( {a-l){a + &) 2 s 
exp er c 



(21) 



To reduce the baryon fraction by factors Cb <; 2 requires 
density variations of magnitude a c > 1. 

In the following discussion, we highlight a set of five 
specific models, listed in Table 1. Models A and B have 
small (~ 20%) baryon corrections, models C and D have 
large baryon bias Cb ~ 2 and model E is an extreme model 
in which the baryon fraction is reduced by an order of mag- 
nitude. 



3.2 Sunyaev-Zel'dovich Effect 

The SZ effect is produced by inverse-Compton scattering 
of cosmic microwave background (CMB) radiation off ther- 
mally excited electrons in the hot ICM plasma (see Birkin- 
shaw 1998 for a recent review). We calculate the central 
Comptonization parameter of the (nonrelativistic) thermal 
SZ effect 



2/(0) = 



k B T(l) 
n £ {l)(jT 5- dl 



(22) 



where the integral dl is along a narrow line of sight through 
the center of the spherical cluster. Here n e — p/ jj, e m p is the 
electron number density and ot the Thomson cross section. 
Since the plasma phases are assumed isobaric, the product 
n e (r)T(r) is constant, and no phase integral is necessary in 
the multiphase case. Deviation of the ^-decrement from that 
of a single-phase plasma is caused by the alteration of the 
overall pressure profile in the cluster. 




Figure 3. The baryon multiphase correction factor evaluated at 
the virial radius, Cf,(r 200) is displayed in a contour plot within 
the e, a c plane. Standard parameters r c /r2oo =0.1 and /3 = 2/3 
are assumed. Labels refer to models whose multi-phase spectra 
are displayed in Figures ^ and ^| 




0.£ 



Figure 4. Contours of the fractional increase in central Sunyaev- 
Zel'dovich decrement for the multiphase models relative to the 
singe-phase case, assuming r c /r2oo =0.1, = 2/3 and a = 0.36. 



The fractional deviation of the central Comptonization 
parameter Ay/y a = [y m (0) — j/ s (0)]/j/ s (0) in the multiphase 
with respect to the single-phase model is shown in Figure ^j. 
In the case of a uniform, multiphase structure (e = 0), the 
fractional change in y follows from the pressure correction 
factor, equation (|l^), 
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A,/, 8 = exp( (1 - Q) 4 (2 - a) ag)-l. (23) 

Figure ^| shows data for the case a = 0.36, approximately 
the slope of the 2—10 keV luminosity versus temperature, 
derived from a Raymond-Smith plasma code assuming a 
one-third solar abundance of metals. Similar values result 
for the case a = 0.5. 

The two models with 20% baryon diminution have mod- 
est, but potentially discernable, SZ effects. Model A has a 
central decrement enhanced by 10% while model B is en- 
hanced by 25% over the single-phase case. The latter is sim- 
ilar to the 30% effect for model C, one of the large baryon 
fraction diminution models. The other factor 2 baryon frac- 
tion model — model D — has a central value of y increased 
by ~ 50% over the standard /3-model . The extreme model 
E has a signal enlarged by a factor 2 over the standard case. 

Even in the single-phase case, there is inherent uncer- 
tainty in predicting the SZ effect amplitude from X-ray ob- 
servations which arises from uncertainty in the physical dis- 
tance to the cluster. At low redshifts, the distance error is 
completely due to uncertainty in the Hubble constant. Given 
a cluster with fixed X-ray properties, a fractional deviation 
in central SZ decrement Ay/y s due to a multiphase medium 
could, instead, be interpreted as a distance effect. This would 
imply a fractional error in the Hubble constant 

AH /H = -2Ay/y s . (24) 

For example, in a universe with true Hubble constant of 
65 km s _1 Mpc -1 , observations of a multiphase model A 
cluster would yield a value 52 km s _1 Mpc -1 and model 
B would produce an estimate of 42. The other highlighted 
models would produce even lower estimates of Ho. 

Note that this result has the opposite sense of correction 
for the SZ decrement compared to other estimates of the SZ 
effect with multiphase gas (e.g.. Holzapfel et al. 1997). This 
is because the gas pressure for isobaric density fluctuations is 
greater than that for a single-phase medium, whereas other 
models with adiabatic density fluctuations, such as those 
present in SPH calculations without cooling (Inagaki, Sugi- 
nohara & Suto 1995), have a pressure lower than that of the 
single-phase gas. 

3.3 X-ray spectra 

Spectroscopic analysis of the X-ray emission provides an al- 
ternative, independent diagnostic of multiphase structure. 
We calculate expected X-ray spectra from the multiphase 
plasma in two ways. We first consider a simple model for 
the X-ray bremsstrahlung continuum from the cluster, us- 
ing an emission function of a purely hydrogen plasma, with 
e(E,T) oc T 2 e ~£t and a Gaunt factor of unity. The 
plasma emission function is then A(T) = KqT 1 / 2 , with Ko an 
arbitrary normalization amplitude. Second, a more detailed 
spectrum is calculated, using the Raymond-Smith plasma 
emission code (Raymond and Smith 1977) with metal abun- 
dances one-third of the solar value (Allen 1973). The former 
approach highlights continuum behavior while the latter al- 
lows the study of the behavior of X-ray lines between 0.5 — 9 
keV. 
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Figure 5. X-ray continuum emission for different multiphase 
models listed in Table 1. A value of o = 0.5 is assumed. All mod- 
els have the same emission weighted temperature of 10 s K in the 
spectral region shown. 

We generate the composite spectrum emitted by the 
multiphase cluster atmosphere by numerically integrating 
the weighted emission over the density distribution at a 
particular radius, then intergrating over the cluster volume 
V — r| 00 . Gas beyond the virial radius r2oo is ignored. 

The appearance of the continuum plasma X-ray emis- 
sion in the multiphase case can differ substantially from that 
of the isothermal, single-phase plasma. In Figure ^, we show 
the 0.05 — 100 keV bolometric X-ray spectra of three mul- 
tiphase models (A, C and E) along with the single-phase 
case. All spectra are normalized to yield the same emission 
weighted temperature of 10 8 K. We assume all phases are 
optically thin. The most important effect on the spectra is 
the appearance of both low-energy (E <C ksTx) and high- 
energy (E 3> ksTx ) enhancements of the spectrum with in- 
creasing magnitude of multiphase effects. This shape change 
arises from the blending of gas at temperatures both below 
and above the fiducial 10 8 K value. In the limit of extreme 
multiphase strength (model E), the bremsstrahlung spec- 
trum approaches power-law like behavior. 

A more complete X-ray spectrum of clusters is code. In 
particular, the use of such a code allows invesigation of line 
emission as a diagnostic of multiphase structure. 

Figure |^ shows the simulated emission, derived from a 
Raymond-Smith code, between 0.1 — 15 keV photon energies 
for a plasma with an assumed metallicity equal to one-third 
of solar abundance. Along with the rise of the low-energy 
continuum, the other prominent effect of increased multi- 
phase structure is the strengthening of low-energy (~ 1 keV) 
emission lines. To highlight line versus continuum effects, we 
plot both zero and one-third solar metallicity predictions for 
the emission for each model shown. 

The complex of lines between 0.5 and 1.5 keV presents 
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E (keV) 

Figure 6. X-ray spectra for different multiphase models derived 
from a Raymond-Smith code. All models have the same emis- 
sion weighted temperature of 10 8 K in the 2 — 10 keV region. 
The luminosity of each model is displaced vertically by arbitrary 
amounts for clarity. Solid lines assume a one— third solar abun- 
dance plasma while dotted lines show the multiphase emission 
assuming zero metallicity. For reference, long dashed lines show 
the X-ray spectrum of the single-phase cluster with no metal 
abundances. 



a useful diagnostic for multiphase structure. Included in this 
region of the spectra are the Fe L-shell lines, as well as H-like 
and He-like emission from N, O, Ne and Mg. For example, 
weak baryon bias models (A and B) are readily distinguished 
by this emission signature, as are models the more strongly 
multiphase models, C and D. 

In contrast to the low energy lines, the strength of the 7 
keV iron complex is almost unaffected by multiphase struc- 
ture. These lines originate in hot phases very close to the 
fiducial temperature of 10 s K. The emission weighted tem- 
perature constraint imposed on the models requires that the 
contribution to the total emission from phases near the fidu- 
cial temperature cannot vary by large factors. Hence the hot 
emission lines do not vary significantly among the multi- 
phase models. 

Given the very different behavior of the low and high 
energy line emission, we investigate the behavior of a hard- 
ness ratio 

L to t[0.6-1.5keV] 



n = 



L tot [6.6-7.5keV] 



(25) 



in the multiphase model plane. For the single phase case, 
1Z = 2.69. Contours of constant 1Z in the {cr c ,e} plane are 
shown in Figure |^. From this figure, it is clear that even 
moderate signal-to-noise spectra could produce useful con- 
straints among models with similar baryon fraction correc- 
tion factors. Models A and B differ in 1Z by nearly a factor 2. 
Models B and C are nearly degenerate in this measure, but 
inspection of Figure ^ shows that model C is more contin- 
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Figure 7. Contours of constant hardness ratio TZ, equation (^) 



uum dominated at low energies while model B has a larger 
line contribution. 



3.4 Limiting the baryon fraction bias 

The concluding sentence of GT96 expresses a view "that 
the intracluster medium is much more complex than most 
people have hitherto assumed and that there is sufficient 
uncertainty in its modeling to permit a critical density, 
Einstein-deSitter universe". Within the context of the ex- 
panded version of their model which we develop here, we 
can ask whether this opinion is supported by recent data. 

A number of high quality measurements of the Hub- 
ble constant from SZ and X-ray observations have been 
made recently (Myers et al. 1997; Birkinshaw & Hughes 
1994; Jones 1995; Hughes & Birkinshaw 1998; Holzapfel et 
al. 1997). Hughes & Birkinshaw (1998) present an ensemble 
value Ho = 47± 7 km s _1 Mpc -1 from these studies. For a 
true value Ho = 65 km s _1 Mpc -1 , supported by Type la 
supernovae (Hamuy et al. 1995; Riess et al. 1996), expand- 
ing photosphere of Type II supernovae (Schmidt et al. 1994) 
gravitational lens time delays (Kundi et al. 1997; Schechter 
et al. 1997), this ensemble average is low by 28%. Assum- 
ing that a multiphase structure is at least partly responsible 
for this biased estimate — other effects may lead to an un- 
derestimate at the ~ 5 — 10% level (Cen 1998) — then a 
bound on the SZ decrement enhancement Ay/y s < 0.14 re- 
sults. Comparing the contours in Figures ^ and [|, this limit 
restricts baryon fraction diminution factors to be modest, 
0.75 < C b < 1. 



3.5 Caveats and extensions 

The model we present contains a number of simplifying as- 
sumptions. It is important to note that the model is, in prin- 
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ciple, falsifiable. Given a known Hubble constant, likelihood 
analyses of SZ observations and X-ray spectra will indepen- 
dently identify preferred regions in the {a c , e} plane. If these 
regions are consistent, the observations can be combined to 
yield a best estimate location {<r c , e}, and an estimate of the 
baryon fraction bias Cb (Figure^). Inconsistent constraints 
may imply a need to relax one or more of the following as- 
sumptions. 

Lack of spherical symmetry. With rare exceptions, cluster 
X-ray images are close to round. Most have axial ratios 
b/a > 0.8 (Mohr et al. 1995). Such small deviations from 
spherical symmetry lead to scatter, but little bias, in deter- 
minations of Ho from SZ+X-ray analysis (Sulkanen, 1999). 
Given supporting evidence for a multiphase ICM in a cluster 
of moderate ellipticity, the spherical model introduced here 
could be extended to prolate or oblate spheroids. A more 
profitable approach might be to include multiphase struc- 
ture in the deprojection method discussed by Zaroubi et al. 
(1998). 

Non-isothermal emission weighted temperature proBles. 
There is indication from ASCA observations (Markevitch 
et al. 1998) that the emission weighted ICM temperature 
declines substantially within the virial radius. However, 
ROSAT colors rule out a temperature drops of 12/20% for 
5/10 keV clusters within one-third of V2oo (Irwin, Bregman 
& Evrard 1999). It is straightforward to include a radial 
temperature gradient Tx{r) into the analysis, entering into 
the definition of the pressure profile, equation 

Non-lognormal distribution of density fluctuations. The cho- 
sen form of the density distribution is motivated by sim- 
plicity and by the observation that non-linear gravity on 
a Gaussian random density field characteristically gener- 
ates a log- normal pdf (Cole & Weinberg 1994). The results 
are sensitive to low order moments of the density distri- 
bution. We await observations and future numerical simu- 
lations including cooling and galaxy-gas interactions in a 
three-dimensional setting to shed light on the appropriate 
form of the density fluctuation spectrum. 

Non-isobaric equation of state. This may be the most read- 
ily broken of our model assumptions. The cluster environ- 
ment is very dynamic. During large mergers, the behavior 
of the gas in the inner regions of infalling subclusters is es- 
sentially adiabatic (Evrard 1990; Navarro, Frenk & White 
1993). During quiescent periods between mergers, a clus- 
ter atmosphere may stabilize and develop the assumed iso- 
baric perturbations during a cooling flow phase (Thomas et 
al. 1986). In the transition period, isobaric perturbations in 
the core may co-exist with adiabatic perturbations near the 
virial radius. Empirical constraints will come from improved 
spectroscopic imaging. 

Binding mass estimates under hydrostatic equilibrium. As 
noted in §2.2, the radial dependence of our model multi- 
phase distribution can lead to a total mass profile M to t,m(r) 
that is steeper than that determined for a single-phase gas. 
Galaxy cluster lensing observations could be used to test 
the mass distribution predicted by multiphase models (see 



Figure |2|) . This provides another independent constraint on 
the admissable region of the multiphase {<j c ,e} plane. 



4 SUMMARY AND DISCUSSION 

We present a spherically symmetric, multiphase model of 
the intracluster medium in galaxy clusters. The model as- 
sumes existence of a lognormal distribution of isobaric den- 
sity and temperature fluctuations at any radius. The radially 
dependent variance of the density fluctuations u 2 (r) is sub- 
ject to two empirical constraints : 1) that the broadband 
X-ray emissivity profile matches observations and 2) that 
the X-ray emission-weighted temperature is constant with 
radius. 

We calculate the bias introduced in cluster gas mass 
fraction estimates when a single-phase model is a applied 
to a multiphase atmosphere. As derived by GT96, the stan- 
dard analysis of the X-ray observations with a single-phase 
assumption will overestimate the baryon fraction in the 
multiphase case. Examining observable effects on the cen- 
tral Sunyaev-Zel'dovich decrement as well as X-ray spec- 
troscopy, we demonstrate how, within the context of this 
model, the bias can be recovered by existing and future ob- 
servations. 

Large values of the clumping factor C, hence large re- 
duction in the cluster baryon fraction are not favored by 
current observations. Models with high values of a c produce 
a nearly power-law X-ray bremsstrahlung continuum and 
bias estimates of the Hubble constant. An ensemble mean 
value 47 ± 7 km s" 1 Mpc -1 (Hughes & Birkinshaw 1998) 
arising from recent SZ+X-ray analysis, when compared to 
an assumed value of 65 km s _1 Mpc -1 , suggests clumping 
only overestimates ICM gas fractions by < 20%. 

Spatially resolved X-ray spectroscopy, particularly of 
line emission in cooler region (0.1-3keV), will provide tests 
of multiphase model. Data from the upcoming AXAF and 
XMM missions will be particularly valuable. 
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